A Model For Late Dark Matter Decay 
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The standard cold dark matter cosmological model, while successful in explaining the observed 
large scale structure of the universe, tends to overpredict structure on small scales. It has been 
proposed this problem may be alleviated in a class of late-decaying dark matter models, in which 
the parent dark matter particle decays to an almost degenerate daughter, plus a relativistic final 
state. We construct explicit particle physics models that realize this goal while obeying observational 
constraints. To achieve this, we introduce a pair of fermionic dark matter candidates and a new 
scalar field, which obey either a Z4 or a U(l) symmetry. Through the spontaneous breaking of these 
symmetries, and coupling of the new fields to standard model particles, we demonstrate that the 
desired decay process may be obtained. We also discuss the dark matter production processes in 
these models. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

There is an abundance of evidence to indicate the ex- 
istence of dark matter (DM) , including its necessary con- 
tribution to both galactic stability and structure forma- 
tion in the early Universe [1-3]. The standard ACDM 
cosmological model, in which cold dark matter (CDM) 
makes up 22% of the universal energy budget, provides 
an excellent description of our Universe. However, little 
is known about the particle properties of dark matter. In 
addition, some problems with CDM are encountered at 
small scales. 

A popular class of CDM candidates is weakly interact- 
ing massive particles (WIMPs). In WIMP models the 
DM couples weakly to standard model (SM) particles, 
which allows for scattering/ annihilation processes. These 
serve to keep the dark sector in thermal equilibrium with 
the visible sector in the early Universe and can, with an 
appropriate choice of coupling, cause the DM to freeze 
out with the correct relic density. Interaction with the 
SM is similarly appealing from a detection standpoint, 
potentially providing both direct [4-6] and indirect [7-9] 
signatures of a given model. Nonobservation of these sig- 
natures allows for constraints to be placed on parameters 
such as particle masses or coupling constants. 

Though the DM mass is unknown, some information 
can be inferred from observations of large scale struc- 
ture. For cold dark matter, structure forms hierarchi- 
cally, with the earliest structures formed on short length 
scales, which can then merge to form larger structures. 
This is to be contrasted with hot dark matter in which 
the largest superclusters form first. Numerical simulation 
has shown the CDM scenario to fit observations well [10], 
while hot dark matter is strongly disfavored. 

The CDM model is not-problem free, however, as it 
tends to overproduce small scale power [10-24]. Simula- 
tions predict cusps in the DM density at the centers of 
galactic halos in conflict with observation. CDM also 
over-predicts the number of dwarf galaxies orbiting a 
Milky Way-sized galaxy by about a factor of 10. Al- 



though simulations do not include visible matter, the 
gravitational potential wells they predict would promote 
a level of star formation not observed. Though these is- 
sues may be partially alleviated by tidal disruption and 
other effects, the small scale power problems of ACDM 
are still poorly understood (see e.g. [25, 26] for recent 
work). Such issues have led many to consider a "warm" 
DM candidate, with a mass of keV scale, intermediate be- 
tween hot and cold dark matter. In this work, as in [27- 
34], we will consider an alternative hypothesis in which 
the usual assumption of a single DM candidate is chal- 
lenged. 

We consider a scenario with two WIMP candidates, in 
which one species is unstable to decay into the other. If 
the mass splitting between the two WIMPS is sufficiently 
small, the decay process will leave the overall halo mass 
unaffected, while giving its constituent DM particles a 
small velocity kick. Such velocity kicks heat the dark 
matter halos and cause them to expand, softening the 
central cusps and disrupting small halos [35-39]. Such 
models are appealing, as they can alleviate the small 
scale structure problems, while retaining the attractive 
features of cold dark matter. 

We shall assume the DM decays predominantly via the 
channel 



X 



x + i- 



(1) 



where \* and X denote the heavier and lighter candidate, 
respectively, and I is some relativistic final state. The 
mass splitting between x* an d X is given by 



Am 



m v «e . 



(2) 



where e <C 1. Abdelqader and Melia [35] have shown the 
dwarf halo problem can be solved for e ~ (5 — 7) x 10~ D 
and a decay lifetime of (1 — 30) Gyr. The work of Pe- 
ter, Moody, and Kamionkowski [36] has demonstrated 
that galaxy cusps can be alleviated for a wider range of 
e and r x * , with the most favored lifetimes in the range 
(0.1 — 100) Gyr. Subsequent work by Peter and Ben- 
son [40] has used properties of galactic subhalos to fur- 
ther constrain the allowed values of e, preferring lower 
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values to those favored in [35]. Dark matter decays may 
be further constrained from analysis of their effect on 
weak lensing of distant galaxies as in [41]. However, at 
present such analyses have only placed limits on models 
with much larger values of e than those considered in this 
work. 

An interesting possibility, from an observational stand- 
point, is a decay mode in which the relativistic final 
state, /, consists of SM particles. This allows the pos- 
sibility of verifying the model, via the detection of parti- 
cles produced by decay in our own Galaxy, or of a diffuse 
flux from decays in halos throughout the Universe. Cur- 
rent astrophysical observations place constraints on the 
allowed parameters, via comparison of the decay fluxes 
with relevant astrophysical backgrounds. Reference [42] 
placed stringent constraints on the decay parameters for 
the case in which I is a photon, while Ref. [43] de- 
rived somewhat weaker constraints for the cases in which 
I = Vv or e ± . For I = 7 or e ± the lifetime is restricted 
to be below about 1 Gyr, while a much larger range of 
lifetimes is permitted for I = Vv . 

The aim of this work is to construct a particle physics 
model which can realize the decaying dark matter sce- 
nario. We shall use the criterion specified by Abdelqader 
and Melia [35] [namely, e ~ (5-7) x 10" 5 and r - (1-30) 
Gyr] as a reference point for these models, but given the 
constraints of Ref. [40] , we will choose the more restric- 
tive value of e ~ 10~ 5 and r ~ (1-10) Gyr. In Sec. II we 
introduce and discuss two possible models for decaying 
dark matter, and outline the DM production mechanism. 
Section III focuses on constraints on the models and the 
available regions of parameter space. We conclude in 
Sec. IV. 



II. A DARK MATTER DECAY MODEL 

In order to construct a model which can achieve this 
decay scenario there are certain criteria that need to be 
satisfied. The first and most important of these is the 
need for two candidates with nearly degenerate masses. 
Second, we need either decay of the parent DM parti- 
cle to light SM final states, or to some new light degree 
of freedom. Third, the process needs to occur on times 
scales relevant for the disruption of structure formation, 
and last, we need some viable DM production mecha- 
nism. WIMP-like scenarios are particularly interesting 
on this front, as WIMPs are populated as thermal relics 
and naturally freeze out in the early Universe with the 
correct relic density. 

Two scenarios will be considered. In the first we im- 
plement a variation of the "exciting dark matter" model 
conceived by Finkbeiner and Weiner [27], which involves 
the addition of a dark sector containing a Dirac fermion 
and a real scalar field to the standard model. The intro- 
duced fields obey a discrete Z 4 symmetry, the breaking 
of which leads to a nondegeneracy of the masses of the 
fcrmion's two Weyl components, and an instability of the 



heavier to decay into the lighter. The scenario considered 
in this work differs from [27] in the values of the model 
parameters chosen; in short, we consider longer decay 
lifetimes. 

As a second example, we consider a generalization of 
the scenario in [27], in which we replace the Z 4 symmetry 
with a global U(l) symmetry, requiring the introduced 
scalar field to be complex. The breaking of this U(\) 
will produce a pseudo-Nambu-Goldstone boson, which 
will serve as our light final state for the dominant decay 
channel. We show that production through interaction 
with the SM is impossible in the U(l) model. The second 
scenario is one example among any number of generaliza- 
tions and extensions to the simple Z 4 model; it is simply 
an illustration that decaying DM can be realized in a 
particle model irrespective of the strength of coupling to 
the SM. 

In Sec. II A we shall explore a model in which SM fi- 
nal states are produced in the dominant decay channel. 
In II B we consider the possibility of completely non-SM 
final states. In II C we discuss production, and in II D 
consider the possibility of \* depopulation. 



A. SM Final States (Z 4 ) 

When searching for a light final state for the process 
in Eq.(l) the obvious place to look is the SM, as the ex- 
istence of particles with masses substantially below the 
CDM mass scale (GeV) is assured. To couple the DM 
to the SM we adopt the model put forward in [27]. Al- 
though this was originally intended as a mechanism for 
explaining the observed INTEGRAL/SPI positron excess 
[9] , with a different choice of parameters the model can 
serve our astrophysical aims quite well. We begin with 
the introduction of a Dirac fermion comprised of the two 
Weyl spinors Xil and X2r, which couple to a real singlet 
scalar <j>. The mass eigenstates for the \i an d X2 fields 
(which we will call x an d X*> respectively) are the DM 
in this model. 

We impose a discrete Z 4 symmetry under which the 
fields transform as 

Xi,2 -> iXi,2 , 

<P^-</>, (3) 

but remain singlets under the symmetries of the SM. This 
allows for the following Lagrangian: 

1 2 
C = -d^4> + J2 xt^Xi - ™dxT l X2r (4) 

i 

- Ai0xT l (xi)r - WtoMZ - V(<j>,H,Hl) + h.c. 

At this stage the two mass eigenstates both have mass 
m D . To lift this degeneracy we need to break the Z 4 
symmetry. We break the symmetry spontaneously down 
to Z2 by allowing <j> to obtain a vacuum expectation value. 
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The Higgs potential is given by 



A 



^\ ,2 



Aft 
4 



(5) 



where H is the SM Higgs doublet. The last term in Eq.(5) 
is included as it is allowed by all the symmetries of the 
theory. Minimizing the above potential with respect to 
both d> arid H we obtain the conditions 



2A (0) 2 - 2/4 + a(h) 2 = 0, 
\ h {h) 2 - 2fi 2 h + 2a{<j)) 2 = 0. 



(6) 



It can now clearly be seen that (0)^0. As in [27], the 
spontaneous breaking of the discrete Z4 symmetry will 
lead to the formation of domain walls, which may be dis- 
favored by observation. We can remove this potentially 
troubling phenomenon by introducing the explicit break- 
ing term /up 3 to our Higgs potential, where /1 is small for 
reasons of technical naturalness. 

Perturbing Eq.(4) about the vacuum there arise Majo- 
rana masses Xi((j>) and \2{4>) for xi an d Xii respectively. 
The Lagrangian therefore contains the mass matrix 



£ D - ( Xu (x 



fer \° \( A l(^) 



2r) 



in,, 
2 



in,, 



(Xil)' 

X2R 



(7) 



which we then diagonalize to obtain the Majorana mass 
eigenstates x ancl X* : 



1 



x ~ 71 + X2)C ~ ^ Xl + X2 ^ 

X* — [(xi + X2) c + (xi + X2)] 



whose masses we find to be 



l x ,x 



= -Jm 2 D +AX 2 _(cb) 2 ± A+(0), 



(8) 



(9) 



where A± = |(Ai ± A|). We want the mass splittings to 
be small, so we choose mn 3> X±(4>), making m x * iX ~ 
^ ± \+(4>), and thus m x -e = 2\ + {(f>). Typically in this 
model we shall consider masses in the range m XiX » ~ (50- 
800) GeV, a breaking scale of (0) - (3-20) MeV, and 
coupling strength of \± ~ 10 _1 [implying Am ~ (0.4-8) 
MeV for e ~ (0.7-1) x 10~ 5 ]. For a detailed discussion of 
the parameter space, see Sec. III. 

In the basis of the mass eigenstates, the Lagrangian 
contains the following interaction terms, which mediate 
both decay and scattering or annihilation processes: 



£ D \+4>xx - A+^x*X* 
- m xXX -m x *~)CX 



A-0X75X* 

- h.c. (10) 



It should be noted that interaction terms coupling like 
mass eigenstates are scalar, while off-diagonal coupling 
is pseudoscalar. As will be seen this has a substantial 
effect on the DM decay rate. 




v, e 



v, e + 



FIG. 1. Primary DM decay channel for the Z4 model 



Mixing of the SM sector with the dark sector (x, X*i 
(j)) occurs through the last term in Eq.(5). Expanding the 
Higgs potential about the vacua of both and H pro- 
duces off-diagonal mass terms for both fields. Expressing 
the potential in terms of the mass eigenstates </>' and h' , 
we find the following mixing of states: 

0~cos6»0'-sin6l/i'~ 4>'-0ti, (11) 
h ~ cos9ti + sin(9 0' ~ ti + 0(j>' , 

where h is the SM Higgs boson, and A^ is its self-coupling. 
To first order in a 



e 



and we find masses of . 



\h(h) ' 
' and h! to be 

-2A (0) 2 , 

— ^h(h) 2 , 



(12) 



(13) 



(in the limit m^ <C m/,) with that of </>' being ~ (2 — 20) 
MeV. In this work we adopt a SM Higgs mass of m^ ~ 
130 GeV. Through 0-Higgs mixing x and x* can couple 
to h! and by extension the SM. In particular, it allows 
the possibility of decay into SM final states via processes 
such as that shown in Fig.l. 

This process has a decay rate given by 



*y?o 2 5 7 

2800 7r 3 mi T7Ve : 



(14) 



where yi is the Yukawa coupling for the dominant SM 
final state, and we assume Am 3> mi, A + ~ A_, and 
that the light final states are Dirac. Just as we choose 
our splitting to be small, we can choose the region of 
parameter space in which the lifetime is sufficiently large 
to disturb structure formation. 

The above decay rate contains several elements which 
naturally lead to suppression and thus a long lifetime. 
First, the decay rate is subject to phase-space suppres- 
sion, as there are 3 bodies in the final state. Second, 
it depends on the Yukawa coupling yi, which, given we 
are only interested in decay into light leptons (e + e~, Vi>), 
will be a small number. It is also dependent on the Higgs 
mixing angle 9, which in turn varies depending on the 
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strength of the coupling a. Although there is some free- 
dom of choice with respect to the value of a, we typically 
take a <~ 10~ 5 , which results in a mixing of 9 <~ 1CP 9 
(see Sec. II C for details). Lastly, pseudo-scalar coupling 
between \ and X* means the decay rate contains a factor 
e 7 (as opposed to e 5 for scalar coupling). The conjunc- 
tion of these factors means that the DM lifetime can be 
long without A± being too small, typically A± ~ 10 _1 . 



B. Dark Decays \U(l)] 

An advantage of SM final states is their detectabil- 
ity. While directly observable consequences are a de- 
sirable model building goal, the nonobservation of the 
signatures of a model can lead to constraints, as will be 
seen in the next section. Were observational constraints 
to strengthen, there is the possibility that a nonobserva- 
tion of the final states in the above model may rule it 
out. Should this occur the viability of decaying DM in 
general would rely on the primary decay channel being 
independent of the SM. In this section we will present a 
model which can realize this. 

One way to naturally produce a decay channel with a 
light final state is to upgrade the discrete Z 4 symmetry 
to a global U(l) symmetry. Spontaneously breaking this 
will produce a Nambu-Goldstone boson (NGB) in the 
theory, which will couple to the DM. The Lagrangian in 
this scenario is similar to Eq.(4) except that now = 
(l/\/2)(0i + i<p2), and \i,2 and <p now transform under 
the [/(!) symmetry as 



Xl,2 



Xi,2e w * , 



(15) 



where 9 X is some arbitrary phase. The Higgs potential 
will have a similar form to Eq.(5) only now <j> ^ <jy , and 
(f> 2 terms become (ftcj). As with the discrete case, we end 
up with the mass matrix in Eq.(7) and subsequent mass 
eigenstates \ and x*; only now they couple to both the 
mass eigenstates 4>'i defined by 



f>i ~ cos 9i 4> x — sin 9\ h! 



a(h) h , 



M0> 



(16) 



h ~ cos9 1 ti + sin9 1 0i ~ ti + 4>[ , 

(m<f, ^> rah as will be seen) and the NGB <j> 2 , i.e. 

It should be noted that the above coupling to the NGB 
is scalar, which results from the fact that fa is the imag- 
inary component of (p. This means that decays into the 
NGB contain less e suppression (one power of e) than 
they would were the coupling pseudoscalar (e 3 ). 

Coupling to the NGB implies the existence of long 
range DM-DM interactions, which can potentially af- 
fect structure formation. To avoid the issues involved 



with this, we will introduce a small soft breaking term, 
^(</> 2 + (ft 2 ) to the Higgs potential to explicitly break 
the continuous U(l) symmetry down to the discrete Z 4 . 
This gives the NGB an O(p) mass, which is naturally 
small. 

The off-diagonal interaction term in Eq.(17) leads to 
the decay channel x* X + 4>2 , which has the decay rate 



47T 



m v *e . 



(18) 



For the case where the primary decay of the DM was 
into SM final states, the reason for a long lifetime was the 
weak mixing with the SM and the suppression from the 
high power of e. Should, however, the decay referred to 
in Eq.(18) be the primary channel, there is no such sup- 
pression, and we are forced to impose the additional ap- 
proximate symmetry Ai ~ — A2 to make A + <~ 0(1O~ 18 ), 
hence achieving r <~ O(Gyr). As Am = 2\ + {4>), small 
A + implies a high breaking scale for reasonable values of 
the DM mass, roughly (<j)} - 10 14 GeV. As will be dis- 
cussed in IIC 2, A0 <~ 1 making ~ (<j>) in this model. 

This high scale has the potential to cause problems. 
Recall the minimization conditions in Eq.(6). In order to 
reproduce the correct breaking scale for the SM Higgs, 
either a needs to be small enough such that a(<f>} 2 is 
negligible with respect to fi 2 , or we have a finely tuned 
scenario resembling the hierarchy problem of the SM, in 
which a((j>} 2 and /i^ are of similar order. The former, 
though the more natural of the two cases, precludes pro- 
duction via mixing with the SM, so we will entertain the 
latter for the time being. 



C. Production 

Both scenarios presented have all the required elements 
to disrupt structure formation in the desired fashion. All 
that is needed now is a production mechanism for the 
dark matter candidate in each scenario. As mentioned 
previously, one of the appealing properties of a WIMP is 
attainment of the correct relic abundance through ther- 
mal freeze-out from the bath in the early Universe. 



1. Z4 Case 

In the model presented in II A, \ couples to the SM 
through the Yukawa sector. It therefore follows that it is 
through these channels that it will maintain equilibrium 
with the SM prior to freeze-out. Production differs from 
the standard WIMP scenario in that it is a two-phase 
process. The <j) is populated via interactions with the 
SM in the Higgs sector, while x and \* are produced 
through their coupling to 4>. At some temperature below 
the DM mass the (j)-x annihilation rate will drop below 
the expansion rate, and \ and x* wm freeze out with 
respect to <fi, fixing the comoving DM abundance to the 
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FIG. 2. Dominant production mechanism for T > run- 



FIG. 3. DM produced through cj> annihilation. 



standard value. We will now chronologically step through 
the processes leading to DM freeze-out. 

The requirement that <f> be in chemical equilibrium 
with the SM well before freeze-out places a constraint 
on the allowed values of the coupling a. Well above 
the electroweak scale, the dominant process keeping cf) 
in chemical equilibrium with the SM will be hh — > <p<f) 
(Fig. 2), which at temperatures well above the Higgs mass 
has the annihilation rate 



T{hh -> (t>4>) 



a 2 T 
256tt 3 



(19) 



This process will remain in equilibrium until T < rrih, 
and h production becomes Boltzmann suppressed, caus- 
ing this cj> production channel to become unavailable. For 
the DM masses of interest in this model (of order or be- 
low rrih), we require a > 10 -6 to ensure that <f> is in 
equilibrium at some point prior to <j)-x freeze-out. 

Also contributing to <j) production is the /i-mediated 
processes // — > <fi(f>, which have the annihilation rate of 



r(// 



a 2 y 2 f W 2 T 3 

167T 3 mt 



(20) 



up to a color factor for processes involving quarks. The 
temperatures at which these processes freeze out depend 
on both the Higgs mixing a and SM fermion Yukawa yf. 
For the values of a considered in this paper (a ~ 1CP 5 ) 
the process which remains in equilibrium longest is that 
involving b quarks. This freezes out around the time 
at which the annihilation in Fig. 2 turns off. Thus the 
temperature at which <f> freezes out with respect to the 
SM can be calculated to be Tf SM ~ 20 GeV. This occurs 
when cf> is still relativistic. 

After <f>-SM freeze out, the temperature of the sys- 
tem will continue to track that of the background. 1 The 



1 Up to a factor (g*/g*) 1//3 , where g* and g' t are measures of the 
number of freedom in the bath at M and ~ x respectively. 
We will assume this ratio to be ~ 1, and y/gZ — 10.8, i.e. that 
all degrees of freedom are in equilibrium. While depending on 
the time of DM freeze out this might not be strictly true, the 
effect on the results will be negligible. It is therefore irrelevant 
exactly when (f> freezes out with respect to the SM, as long as it 
has been in equilibrium at some point prior to freeze out. 



DM will be kept in chemical equilibrium with tf> through 
the scattering in Fig. 3, which in the nonrelativistic limit 
has a cross section of 



av re l 



|A+| 4 



(21) 



This process will freeze out once the temperature of the 
4>~X system falls below m x and the number density of x 
becomes Boltzmann suppressed. 

To determine the DM relic abundance we use the well 
established result [441 



o aXDM-f \fqZGeV 1 
n^h 2 = 1.07 x 10 9 /yy , ; 

g* s mpi{av) 



(22) 



where 



XDM-f =m x /Tj x = In[0.038(# /y/gZ)mpim x (<Tv)] 

--]n[ln[0m8(g/^)mpim x (av}}} , (23) 

and Tf~ x is the temperature at which cf> and x drop 
out of chemical equilibrium. We find that typically 
XDM-f ~ 20. The requirement that we produce the ob- 
served relic density of ft x ~ 0.22, places constraint on the 
free parameters A+ and m x should the DM be a thermal 
relic. See Sec. Ill for a full treatment of the parameter 
space. Given that m x ~ m x * (and Am -c Tp x ) x an d 
X* will be produced in equal abundance. 

After <j)-x freeze-out, the relativistic 4> will remain with 
fixed abundance until the spontaneous breaking of the Z4 
symmetry (at MeV scale). After symmetry breaking they 
become unstable to decay into photons via the loop order 
process depicted in Fig. 4 [45, 46]. This process has a rate 
of 



-> 77) 



(24) 



4.5 x lOV 1 



10- 



(10MeV\ 

V m 4> ) 



which is large compared to the expansion rate, and <p 
rapidly depopulates. 

In our calculation of the process depicted at tree level 
in Fig. 3, we have omitted the contribution from lad- 
der diagrams involving (f> exchange in the initial state. 
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FIG. 4. 1-loop order decay cj> —¥ 77, through h-cj> mixing. 
Includes contribution from loops involving W ± , unphysical 
charged Higgs components h , and Fadeev-Popov ghosts. 



This approximation is valid at high energies, but begins 
to break down near freeze out, when the DM is in the 
moderate-nonrelativistic regime. 

At low velocity the Yukawa potential (resulting from 
<f> exchange) from one initial state x can significantly dis- 
tort the wave-function of the other from that of a free 
particle. This leads to an enhancement of the velocity 
averaged cross section in an effect known as Sommerfeld 
enhancement [28, 47-49]. This effect can be taken into 
account by multiplying the relevant cross section by a 
velocity dependent Sommerfeld factor S. To calculate 
the enhancement to the process in Fig. 3 we follow the 
method of [47, 48], but find that in the relevant region 
of parameter space S is close to 1 and the enhancement 
negligible. The enhancement generally becomes more im- 
portant for larger values of Am. 



2. U(l) Case 

Production in the second model presented is slightly 
more difficult. An unfortunate consequence of a small 
Yukawa coupling is a weakening of the annihilation cross 
section (Fig. 3). This suppression ensures that the pro- 
cess in Fig. 3 is never in equilibrium, making thermal 
production of the DM impossible. This leads us to con- 
sider a nonthermal production mechanism, in which x 
and x* are produced out of equilibrium through their 
weak mixing with the bath. Another possibility is pro- 
duction through direct coupling of the DM to the infla- 
ton. While this is clean in that it is independent of SM 
processes, it requires fine tuning to attain the correct 
relic abundance. For the time being we will entertain the 
former possibility. 

The dominant channel through which production can 
occur is through the SM Higgs annihilation pictured in 
Fig. 5. As (4>) is large in this model, we expect this 
process to be strongest above the electroweak breaking 
scale. At these high temperatures finite temperature ef- 
fects come into the Higgs potential at loop order [50-52]. 
This has the effect of giving the scalar components of the 
SM Higgs doublet temperature-dependent masses of the 



FIG. 5. Dominant DM production mechanism in the U(\) 
model. 



form 



KT 2 
24 



(25) 



The process in Fig. 5 goes to a maximum near the (f>i res- 



onance, in which region mi 



Am 2 . 



Granted a <C \ 



and following the analysis of [53], the velocity averaged 
cross section can in this region be well approximated by 



(av) 



(26) 



(^^Jm 2 -Am 2 + A m ) coth (^) + 128A^m 

To avoid Boltzmann suppression in Eq.( 26) we will 
take A<£ to be small for now (A^ ~ 10~ 17 ), which implies 
that m ( p 1 is far below the U(l) breaking scale (m < p 1 ^100 
TeV). In order to calculate the abundance at a particu- 
lar temperature, we must solve the comoving Boltzmann 
equation, which can be expressed in the form 



dn x (T) _ 
dT 

and has the solution 



0\2 



K) 

HT 



{(TV 



n x {T) = T 3 



K) 2 M 

HT' 



dT' , 



(27) 



(28) 



where T n —q ~ \/6/\h is defined by the temperature at 
which m\ x ~ an d will be taken to be when signifi- 
cant production starts. 

For the representative region of parameter space A^ ~ 
1CT 17 , a ~ 1(T 15 , A+ - 1(T 18 , and r x - ~ 1 Gyr, the co- 
moving number density can be calculated to be O(10 -27 ), 
roughly 17 orders of magnitude below the required value 
at that temperature. These values for the parameters 
in the model were chosen as they were shown to maxi- 
mize production. As this channel is expected to be the 
strongest available it is therefore clear that production of 
the DM via mixing with the SM in such a model is im- 
possible. The implication of neither a SM final state nor 
SM related production is the independence of the dark 
sector from the visible. This gives us complete freedom 
in the choice of dimensionless parameters a and A,* but 
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precludes entirely the possibility of direct verification of 
the model. We can now choose ~ 1 and a to be very 
small to avoid issues of fine-tuning. 

Independence of the dark sector from the SM implies 
the necessity for some novel DM production mechanism. 
As mentioned earlier this can be realized through a direct 
coupling of the DM to the inflaton, but as stated such 
a mechanism is problematic as it is difficult to obtain a 
relic density of order that of the SM without fine-tuning 
of the DM-inflaton coupling. 

D. Depopulation of the Excited State 

In the Z4 model, as the temperature of the 4>-\ system 
drops well below the DM mass, x an d X* will have chemi- 
cally frozen out fixing the relic abundance. The s-channel 
equivalent to Fig. 3 will, however, maintain kinetic equi- 
librium in the (f>-x system to temperatures down as low 
as m^ [47, 54]. Both \ an d X* will be kept in equilib- 
rium with each other by way of a process like that in 
Fig. 6, causing both to track closely the temperature of 
the background. However, as the average kinetic energy 
drops below Am the process XX ~~ ^ X*X is n0 longer 
kincmatically viable, and the up-scattcring rate becomes 
Boltzmann suppressed [54]. The result is a rapid de- 
population of x*, and an absence of the heavy state so 
important for disturbance of structure formation. This 
issue can be averted should the scattering rate for Fig. 6 
be small enough such that the process freezes out suffi- 
ciently early, i.e for T ^> Am. Should this be the case, 
both the forward and back scattering processes will cease 
well before depopulation becomes an issue. 

The cross section for this process (at tree level) can be 
calculated to be 

3|A_A+| 2 I [ 32 1 

7rm£, V re l L<J 

in the limit m^ <C m x »Am, which is justified in the re- 
gion of parameter space considered (see Sec. III). In the 
moderate to nonrelativistic regime, the scattering rate 
for process x*X ~ ► XX is given by 

T ~ (n x )^jL [\av rel )Sv 2 rel e- x >*-f<'/ 4 dv rel (30) 
V 71 " Jo 

where x sc .f — m x /Tj : ~ x , and Tj~ x is defined as the 
temperature at which the process in Fig. 6 freezes out. 

After the process in Fig. 3 freezes out, the comoving 
DM number density n x /T 3 is fixed, and is given by 

n x /T 3 ~ ^ 3.76 x ItT 11 GeV . 

m x 

where T is the temperature of the bath. We can now 
choose parameters such that the process in Fig. 6 freezes 
out around the same time as that of Fig. 3, in which 
case Eq.(21) and Eq.(29) are of a similar order. In the 
relevant region of parameter space, this is generally the 




FIG. 6. Process by which \ an d X* maintain chemical equi- 
librium. 



case, with x sc _f ~ 1. Interestingly, this is before tfi-x 
freeze-out, meaning x an d X* are both are in equilibrium 
with (f> but not each other. 

In the above, we have considered only the depopulation 
of x* i n the early Universe. It is also important that x* 
not be depopulated via scattering in the late Universe, 
when Sommerfeld effects are significant. In fact, addi- 
tional constraints on XX an d XX* scattering arise from 
the requirement that self-scattering of DM does not sig- 
nificantly perturb galactic halo shapes [55]. These re- 
quirements will be taken into account in Sec. III. 

In the U(l) model there are no such depopulation is- 
sues, as A+ is very small and the DM is never in equilib- 
rium in the first place. 



III. CONSTRAINTS ON Z 4 MODEL 

Up to this point there has been minimal discussion of 
the choice of values for the many free parameters in our 
model. In order to do so clearly it is important to un- 
derstand exactly what constraints are present. There are 
initially 7 independent free parameters, those related to 
the fermions x an d X*i namely, m x * and A±, and those 
belonging to the Higgs sector: A^, a, /j,^, and [i. Recall 
also that we can express the mass splitting in terms of 
these parameters, that is, Am = 2\ + ((f>) [((f)) depends on 
Higgs potential parameters from Eq.(6)]. Thus when we 
parametrize Am in terms of e (Am = m x *e) and fix its 
value to e = ICP J , we place a constraining relationship 
between A+, m x », and ((f)). As a second constraint we 
will impose A + ~ A_, as they will only differ greatly in 
the finely tuned scenario where Ai ~ ±A2 to high preci- 
sion. Last, we must satisfy the condition in Eq.(22), to 
ensure correct relic abundance. These three constraints 
reduce the number of free independent parameters to 4, 
which can be taken to be (<fi), \<p, a, and /x. We can now 
express allowed values of Am as a function of breaking 
scale (</>) for chosen vales of a and A^. We must choose a 
appropriately such that <f> goes into equilibrium with the 
bath before the temperature of DM frcczc-out. The al- 
lowed values of Am for the appropriate DM lifetimes are 
plotted on the left hand side of Figs. 7- 8, while the corre- 
sponding values of m x » (for e = I0~ 5 and e = 0.7 x I0~ D 
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for Figs. 7 and 8, respectively) are plotted on the right. 

The presence of readily detectable charged particles 
in the final state increases the possibility of both direct 
and indirect detection. Indeed heavy constraints can be 
placed on the parameter space based on nonobservation 
of the consequences of such a final state. In [42, 43] 
detailed analyses of the photon, positron and neutrino 
backgrounds were performed with decaying DM models 
in mind, and constraints placed on the relevant param- 
eters tdm and Am. We have translated the constraints 
on decay to e + e~ to the parameter space relevant to this 
model, resulting in the exclusion region in Fig. 7. 

Which leptons will be produced predominantly will de- 
pend not only on the choice of parameters (i.e. for which 
lepton docs Am > 2m; hold) but also the choice of neu- 
trino model. We will consider three distinct possible final 
states: (i) e + e~, i.e. Am > 2m e for Dirac neutrinos, (ii) 
vV, i.e. Am < 2m e for Dirac neutrinos, and (iii) vv for 
Majorana neutrinos. 

(i) If we consider the SM neutrino to be a Dirac par- 
ticle, then the upper bound on light neutrino masses 
implies a Yukawa coupling of y v < 1CP 11 [56]. Thus 
when decays into charged leptons are kinematically al- 
lowed (Am > 2m e ), their relatively large coupling in the 
Yukawa sector will render decays into neutrinos subdomi- 
nant. There is the important constraint that Am < 2m M , 
as should pairs be produced, their Yukawa is large 
enough that for no allowed values of Am and (4>) would 
t x » >0.1 Gyr. Thus for Am > 2m e , decays to e + e _ will 
dominate. A representative region of parameter space 
can be seen in Fig. 7. To obtain the correct relic abun- 
dance, parameters must lie on the dashed line. We find 
that for m x ^600 GeV, the breaking scale (4>) is required 
to be in the ~10 MeV range, while (for e = 10 -5 ) Am 
is in the MeV. It should be noted that these parame- 
ters coincide with an x sc .f ~ 1, which is well above 

T?~ x ~ Am, removing the possibility of depopulation 
of the heavier DM state. 

Interestingly, the (1-30) Gyr lifetime range preferred 
by Abdelqader and Melia [35] has been nearly completely 
excluded for decays into charged particles, leaving only 
the restrictive region of (0.1-1) Gyr available. It should 
be noted however that should decays to neutrinos domi- 
nate, we can avoid this exclusion region entirely. 

(ii) Should Am < 2m e only neutrinos are kinemat- 
ically available. As Dirac neutrinos couple only very 
weakly with the SM Higgs, the lifetime of the DM will 
be too long to affect structure formation. There are two 
ways in which we could reduce t x * , by either increasing a 
or decreasing m^. We find, however, that for m^ > Am, 
there are no values of a and m^ that can yield a life- 
time short enough. If, however, m^ < Am, the process 
X* — > X + 4>' becomes kinematically allowed. The rate 
for this process does not contain the high level of sup- 
pression that decays into SM final states suffer, and we 
find its lifetime to be 0.1 Gyr, dominating over decays 
into vv. Thus for the choice of parameters m^ > Am 
the DM lifetime is too long, and for m^ < Am vv fi- 




<0> [GeV] 

FIG. 7. Available parameter space for decays into e + e~ (yi = 
2/ e ) for = 1 and a = 10 -5 for lifetimes r x * = 0.1 Gyr (solid 
black upper line), and t x * — 1 Gyr (solid black lower line). 
Parameters yielding correct freeze-out abundance lie on the 
dashed black line. Shaded is the exclusion region from [43]. 
We have chosen e = 10 . 



nal states are unimportant, and r x * is far too short. It 
therefore seems that in no region of parameter space can 
decays into Dirac neutrinos affect structure formation. 

(iii) Should we introduce Majorana masses for the v R 
and employ the type I see-saw mechanism, we have the 
freedom to make y v large enough (while still keeping the 
neutrino mass small) such that decays with neutrino fi- 
nal states will dominate without the need for fine-tuning 
m-0. We can consider three options: y% <§; y\, y% ~ y\, 
and yl ^> y 2 e . Should yi ~ y\ or yi <§; y\, decays into 
electrons are either important or dominate, and so al- 
lowed parameters will be the same as for the Dirac case. 
However, for ^> yl, neutrino final states are dominant 
for all values of Am, and while we still need to respect 
the observational constraints in Fig. 7, we have a wider 
parameter space available, an example of which can be 
seen in Fig. 8. 

Conversely to before, having Am > 2m e makes r x * > 
0.1 Gyr impossible, as the Yukawa controlling the decay 
is much larger than that of the electron. This constraint 
requires us to choose much smaller values of than for 
e + e~ final states, if we wish to maintain thermal pro- 
duction. The smaller the value of m^, the larger the x~X 
cross section in present-day halos. In [55] authors argued 
that to maintain the observed ellipticity in galactic ha- 
los, the timescale for DM self-interactions must be longer 
than the halo age (Pdm-dm > ^ W ^ r )- Following the 
approach in [55] constraints were placed on our param- 
eter space, resulting in the shaded exclusion region in 
Fig. 8. 



IV. CONCLUSIONS 

Models for decaying dark matter are interesting in 
that they maintain the attractive features of the ACDM 
model, while alleviating the issues pertaining to the over 
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<0> [GeV] 

FIG. 8. Available parameter space for decays into ~Dv with a 
larger Yukawa of y v ~ 10~ 4 (y are Majorana) for \^ = 0.8 
and a = 1CP 5 for lifetimes r x « =0.1 Gyr (solid black upper 
line), t x » = 1 Gyr (solid black center line) and t x * = 10 Gyr 
(solid black lower line). Parameters yielding correct freeze- 
out abundance lie on the dashed black line. Shaded is the 
exclusion region based on ellipticity constraints [55]. We have 
chosen e = 0.7 x 10" 5 . 

prediction of small scale power. In this work we investi- 
gated two examples of the class of DM models in which 
decay occurs via the process x* X + U where x* an d 
X are nearly degenerate in mass (in this work we chose 
Am/m x * = e ~ 10 -5 ) and I is relativistic. 

In the first scenario, we considered the possibility of de- 
cays into SM final states. We demonstrated that through 
the breaking of a discrete Z4 symmetry with the real 
scalar field <j>, we could both produce two Majorana DM 
candidates x* an d X with nondegenerate mass, and al- 
low for the decay channel x* X + SM. The required 
long lifetime [(0.1-100) Gyr] was naturally achieved, as 
the the decay rate was suppressed by a high power of 
e, by small Yukawa couplings, and by the small mixing 
between SM-sector and dark-sector particles. 

The only two viable decay modes involving SM final 
states were x* ~ > X + e + e~ and x* ~~ ^ X + v y i where 
the latter is possible only in the case of Majorana neutri- 



nos. We found that for DM masses in the range (50-800) 
GeV [Am ~ (0.4 - 8) MeV] and for a 0-Higgs mixing 
of a ~ 10~ 5 , all required criteria, including thermal pro- 
duction, could be met if the Z4 symmetry was broken at 
the MeV scale, with ((/>) ~ (3 - 20) MeV (m ~ (2 - 20) 
MeV). Interestingly, in applying the constraints on de- 
cays to e + e~ derived in [43], we showed that this final 
state is almost excluded for DM lifetimes in the (1-30) 
Gyr range preferred in [35]. Dirac neutrinos were un- 
able to fulfill the requirements for decaying DM, as their 
Yukawa coupling is too small. Thus decays to Majorana 
neutrinos are preferred by such DM decay models, as 
they are not constrained to either the short lifetimes ap- 
plicable for e + e~ decays nor the small Yukawa couplings 
of Dirac neutrinos. 

In the second scenario, we considered the possibility 
of non-SM final states. This was achieved by replacing 
the discrete Z4 symmetry with a continuous U(l) sym- 
metry. Breaking of the U(l) symmetry led to a pseudo- 
Nambu-Goldstone boson, which became the light final 
state produced in decays. As the DM decay process was 
no longer strongly suppressed, we were forced to finely 
tune model parameters to obtain a DM lifetime in the 
correct range. A consequence of this fine-tuning was to 
make DM production via mixing with the SM no longer 
possible. In this scenario, the dark and visible sectors are 
almost decoupled from each other. Though aesthetically 
less appealing, this model demonstrates the feasibility 
of decaying dark matter, independent of the strength of 
coupling to the SM. 
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